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The nonholonomic constrained system with second-class constraints is in- 
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vestigated using the Hamilton-Jacobi (HJ) quantization scheme to yield the 
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complete equations of motion of the system. Although the integrability con- 




ditions in the HJ scheme are equivalent to the involutive relations for the first- 



class constrained system in the improved Dirac quantization method (DQM), 



one should elaborate the HJ scheme by using the improved DQM in order to 



obtain the first-class Hamiltonian and the corresponding effective Lagrangian 



having the BRST invariant nonholonomic constrained system. 
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I. INTRODUCTION 



Recently, several interesting constrained systems were investigated |1J in the framework of 
the improved Dirac quantization method (DQM) Based on the Caratheodory equivalent 
Lagrangians method ||, an alternative Hamilton- Jacobi (HJ) quantization scheme for the 
constrained systems was also proposed Q and this HJ scheme was exploited to quantize 
singular systems with higher order Lagrangians such as the systems with elements of the 
Berezin algebra |JJ and the Proca model ||. One of the most interesting application of the 
HJ quantization scheme is the systems with second-class constraints [0,§J, simply because 
the corresponding set of equations is not integrable || . If the singular system is transformed 
to become completely integrable, the Hamiltonian has the form suitable for application of 
the Hamilton- Jacobi equations. 

In this paper we improve the HJ quantization scheme for the nonholonomic constrained 
system (NHCS) studied in the literature |9|,|T0f to obtain the complete solutions of the HJ 



partial differential equations (PDEs) for the system and to compare them with those of the 
standard and improved DQMs. In section 2 we briefly recapitulate the HJ quantization 
scheme. In section 3, in this refined HJ scheme, the NHCS with nonholonomic primary 
constraint is reanalyzed to yield the complete solutions. In section 4 we treat the NHCS 
by using the standard DQM and in section 5 we construct the first-class Hamiltonian and 
the first-class effective Lagrangian corresponding to the integrable system in the improved 
DQM. Moreover, with the first-class effective Lagrangian, we have constructed the BRST 
invariant NHCS. 

II. HAMILTON- JACOBI QUANTIZATION SCHEME 

In this section, we briefly recapitulate the HJ quantization scheme @,[|. We start with 
an unconstrained system with the Lagrangian L, for which we can obtain a completely 
equivalent Lagrangian described as 
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L' = L(^-^§4 (2.1) 

with i = 1, 2, n. These Lagrangians are equivalent to each other if there exists a function 
S(qi,t) such that the Lagrangians L and L' have an extreme value of the action simultane- 
ously. To guarantee this equivalence, one needs to find functions ai(qj,t) and S(qi,t) such 
that, for all neighborhood of <ji = ai(qj,t), 

L'(q i ,q i = a i (q j ,t),t) = 0, (2.2) 

and L'(qi,qi) is positive to yield at % = a<i(qj,t), 

f = 0. (2.3) 

Note that the Lagrangian L' has now a minimum at = ai(qj, t) so that the solutions of the 
differential equations given by % = oti(qj,t) can yield the extremal action. Now exploiting 
Eqs. fl2.1|) and ( f2~2"D one can obtain at = on 

dS T dS 

w- L ~w, q - (2A) 



Similarly, combining Eq. (|2.1| ) and (|2.3|) yields at qi = a.i the HJ equation 

dS , , 

7^-=P„ 2.5 

where are the conjugate momenta. Inserting pi into Eq. ( j2.4| ), one can obtain the HJ 
PDE in terms of the Hamiltonian H as follows 

dS 

— = -H = -pitji + L(qi, q t ). (2.6) 

Next, we consider a constrained system in which canonical variables are not all indepen- 
dent. In the constrained system, the Lagrangian L is singular so that the determinant of 
the Hessian matrix = Q ^g^. is zero and the accelerations of some variables qi are not 
uniquely determined by the positions and the velocities at a given time. 

Now we consider the rank n — m of the Hessian where determinant of a sub-matrix of 
the Hessian is not zero and thus some velocities q a (a = 1, 2, ...,n — m) can be solved as a 



function of coordinates and momenta p a to yield q a = q a (QitPb)- The remaining momenta 
p a (a — n — m+ 1, n) are functions of qi and p a to yield 

p a = -H a (qi,p a ), (2.7) 

which are equivalent to the primary constraints p a + H a in the Dirac terminology . The 
Hamiltonian Q2.tj|) then becomes 



H$ = p a q a + p a q a - L(qi,q a ,q a ), (2.8) 

which can be shown not to depend explicitly on the velocities q a . 

With the redefinition: tg_ = (to, t a ) = (£, g a ), (a = 0, n — m + 1, n) and po = |f > Eqs. 



]6|) and ( |2.7| ) yield the generalized HJ PDEs for a — 0, n — m + 1, n, 

H , & = p* + H SL (tp,q a ,p a ) = 0. (2.9) 

Exploiting Eqs. ( |2.8| ) and (p.9|), one can obtain 

rf<3i = -7r^dt a , dpi = — 7r=dt a . (2.10) 
<% - - % - 

Here note that we have used the extended index i [i = 0, 1, ...,n), instead of the index a 
(a = 1, 2, n — m) used in the literature ||, to obtain the complete solutions to the system. 
Eq. (FTIOD then yields 

d5= (-^a+Pa^J^, (2.11) 

from which we obtain the action of the form 

S = J(-H dt + Pi d qi ). (2.12) 

Note that at the moment dqi cannot be integrable to yield the desired effective Lagrangian, 
which will be realized by introducing auxiliary fields in the improved DQM in the next 
sections. 

Now, in order to discuss the integrability conditions, one can introduce a linear operator 
Xg_ {ql — 0, n — m + 1, ...,n) corresponding to Eq. ( [2.10|) as 



Xa - f - {LH - } ~di-dW-W^ (2 - 13) 

from which one can obtain the bracket relations among the linear operators Xg_ 

[Xp,X^f = {f,{H'„,H'p}}. (2.14) 

Note that, if one can introduce operators X a with an extended index a (a = 0,n — m + 
1, ...,n, ...) such that these operators satisfy a closed Lie algebra 

[X s ,X a ]f = {f,{H' &1 H' & }} = 0, (2.15) 

then the system of PDEs Xpf = is complete and the total differential equations dqi = 
fifidtp is called integrable. Since the total differential for any function F can be written as 
dF = {F, H'p}dtp, the integrability conditions for a — 0, n — m + 1, n, are given as 

K = {HI H' } + {HI H'p}qp = 0. (2.16) 

Note that the definition of the brackets (whose index i runs from to n) in Eq. ( |2.16| ) 
slightly differs from that of usual Poisson brackets (whose index i runs from 1 to n). If 
Hp does not possess time-dependence explicitly, the integrability conditions ( 2.16 ) are then 



equivalent to the consistency conditions in the DQM and the involution relations in the 
improved DQM, which will be discussed in the next sections. 

III. NHCS IN HAMILTON- J ACOBI SCHEME 

In this section, we consider the nonholonomic constrained system (NHCS), where the 
primary constraint cannot be expressed in terms of the coordinates only, by introducing the 
Lagrangian of the form P,|i~0| 

with the canonical momenta 

Pi = 4i, P2 = ^(93-^2) + qi + q 3 , p 3 = hq 2 - 4s)- (3.2) 



Since the rank of the Hessian matrix Hij = 1,2,3) is two, we have two independent 
relations of the momenta p\ and p 3 , and the dependent one p 2 given as 

P2 = ~P3 + qi + Q3 = ~H 2 , (3.3) 

which is a nonholonomic primary constraint in the Dirac terminology . 
On the other hand, the Hamiltonian given as 

H = l -(p\ - 2pj) + qi + q 2 + ql (3.4) 

and Eqs. Q and Q yield the generalized HJ PDEs for H'^ (a = 0, 2) 

H' Q = Po + H = 0, H' 2 = p 2 + P3 - qi - qz = 0. (3.5) 



Since the Hamilton equations are given by Eq. ( |2.10| ), the above H' a (a = 0,2) functions 
generate the following set of equations of motion 



dqo = dt, dqi = p\dt, dq 2 = dq 2 , dq 3 = —2p 3 dt + dq 2 , 

dp = 0, dpi = —dt + dq 2 , dp 2 = —dt, dp 3 = —2q 3 dt + dq 2 . 

Note that, since dq 2 is trivial, one could not obtain any information at this level. 
Next, for the above H and H' 2 , the integrability conditions ( |2.16| ) then imply 



(3.6) 



H' = {H' , H' 2 }q 2 = -H' 3 q 2 = 0, H' 2 = {H' 2 , H' } = H' 3 = 0, (3.7) 
with H' z being a nonholonomic secondary constraint of the form 

H f 3 = 2p 3 - Pl -2q 3 -l. (3.8) 
This H 3 then yields an additional integrability condition 

H' 3 = {H 3 , H' Q } + {H' 3 , H' 2 }q 2 = 4p 3 -4q 3 + l-q 2 = 0, (3.9) 
to arrive at the desired information of dq 2 absent in Eq. (|3.6|). 



q 2 -2q 2 + 2 = 0, (3.10) 
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so that we can now solve the equations of motion completely. Here one notes that using Eq. 
]2]) the nonholonomic constraint (P-8Q can be rewritten in terms of the and as below 



H' 3 = -q 1 + q 2 -q 3 -2q 3 -l. (3.11) 



With the aid of Eq. ( |3.10|) , we can now completely find the solutions for the equations 
of motion (R.6I) as 



qi (t) = Ae 2t -t + Cx, Pl (t) = 2Ae 2t - 1, 

q 2 (t) = 2Ae 2t + t + C 2 , p 2 (t) = -t + C 1 , (3.12) 
q 3 (t) = je 2t + Be- 2t + §, p 3 {t) = \Ae 2t + Be' 21 + §, 

where A, B, C\ and C 2 are arbitrary constants of integration. Note that the results (|3.12|) 
are exactly same as those of Ref. @.|1(J except the existence of C\ in our solutions. 

IV. NHCS IN STANDARD DIRAC QUANTIZATION METHOD 

In this section, we analyze the Hamiltonian structure of the Lagrangian (|3.1|) in the 
standard DQM to compare with that in the HJ scheme. With the definition of the canonical 
momenta (|3.2|) one can obtain the nonholonomic primary constraint of the form, which is 
the same as H' 2 defined in Eq. (|3] 



^i = P2 + Ps ~ qi - Qs ~ 0. (4.1) 
Now we define the Hamiltonian H with a Lagrangian multiplier v, 

H = H + vn 1 , (4.2) 



from which, requiring the time stability of the primary constraint ( |4.1|) , one can easily find 
secondary constraint, which is equal to H' 3 in Eq. (|3.8| ), 

= {fix, H} = n 2 = 2p 3 - Pl - 2q 3 - 1 w 0. (4.3) 

The time stability of fl 2 yields 



h 2 = {n 2 , H} = 4p 3 -4q 3 + l-v = 0, (4.4) 

to fix the v as v = 4p 3 — 4g 3 + 1. These constraints in Eqs. (4.1) and ( f4.3|) make the system 
second-class with A ab = {Q a , Qf,} = e a b and E\ 2 = 1- 

Next, to obtain the equations of motion for q^ (i — 1, 2, 3), we can proceed to construct 
the Poisson brackets {qi, H} and {pi, H} for the physical variables (qi,Pi) to yield 
qi = qi, q 2 = 4p 3 - 4q 3 + 1, q 3 = 4p 2 + 6p 3 - 4qx - 8q 3 + 1, 

(4.5) 

Pi = Ap 3 - Aq 3 , p 2 = -1, p 3 = Ap 2 + 8p 3 - 4gi - 10q 3 + 1, 

which, together with the constraints (|4.1|) and (|4.3|) , reproduce the solutions (|3.12|) . Note 
that the equation of motion for q 2 in Eq. (|4.5|) is the same as the above fixed value of v. 
Note that the Poisson brackets in the HJ scheme are the same as those in the DQM since fl a 
do not depend on time explicitly. Moreover, if the integrability conditions ( |3.7|) and ( p.9|) 
are rewritten in terms of f)i(= H' 2 ) and fl 2 (= H' 3 ) and v(= q 2 ), one can easily reproduce 
Eqs. ( |4.3j ) and ( |4.4| ) to explicitly show that the integrability conditions in HJ scheme is 



equivalent to the consistency conditions in DQM. 

On the other hand, in order to consistently quantize the NHCS, one should obtain the 
Dirac brackets 

{Qi,Pj}D = 5ij - 5a((5ji + 5 j3 ) - 25 i2 S j3 + 25 i3 8 j i ) 

{Qi, Qj}d = -2(5 i2 8j 3 - 5 i3 5 j2 ) - (8 il 5j 2 - S i2 8j 1 ) + (8 i3 8ji - 8n8 j3 ), 

{Pi,Pj}D = -2(8n8 j3 - SaSjt), (4.6) 

where the Dirac brackets for any functions A(q, p), B(q, p) are defined as {A, B}d = {A, B} — 
{A, Q a }A ab {Q b , B}, with A ab being the inverse of A ab . 

V. NHCS IN IMPROVED DIRAC QUANTIZATION METHOD 
A. Gauge invariant Lagrangian for NHCS 



Now, according to the improved DQM we embed the second-class constrained 

system into first-class one via systematic first-class prescription, where one introduces an 

8 



auxiliary canonical pairs (9, Tie) satisfying {9, tcq} = 1 to yield modified first-class constraints 
Q a satisfying a closed Lie algebra {Q a ,Cl b } = 0. Following the improved DQM |I],|2| we can 
find effective first-class constraints as 

= + 9, n 2 = n 2 -7i e , (5.1) 

whose Poisson brackets strongly vanish in the extended phase space due to the introduction 
of the auxiliary canonical pairs of (6,110). Note that, in the limit (9,ttq) — > 0, Q\ and fl 2 
reduce into H 2 and H 3 in Eqs. Q3.5p and ( |3.8| ) of the HJ scheme, respectively. On the other 
hand, we can also obtain first-class physical variable in this first-class embedded phase space 
as 

q 1 = q 1 - Q q 2 = q 2 + n e , q 3 = q 3 + n e + 29, 

(5.2) 

Pi = Pi + 7r e , p 2 = p 2 , p 3 = p 3 + Tie + 29. 

Note that the Dirac algebra ( (4.6|) in the original phase space is mapped into the Poisson 
algebra of the first-class physical variables in the extended phase space. Moreover, since 
the second-class nature of the NHCS is sometimes suffering from unfavorable problems such 
as ordering upon quantization, it is preferred to convert the second-class system into the 
first-class one so that one can perform consistent quantization. Using the first-class physical 
variables ( |5.2| ), we can now construct the first-class Hamiltonian as 

H = H + (-4p 3 + 4g 3 - 1)9 + (pi ~ 2p 3 + 2q 3 + 1)tx 6 + \ 2 d , (5.3) 

where Hq is the original Hamiltonian (|3.4|) . Moreover, one can also construct the equivalent 
first-class Hamiltonian of the form 

H' = H + 7r e Q 2 , (5.4) 

to satisfy the Gauss law constraints {Qi,H'} = Q 2 and {f^)-^'} = 0. Note that these 
constraint Lie algebra are the same as those in the HJ scheme in the (9, tyq) — > limit to 
show that the involution relations are equivalent to the integrability conditions in the HJ 
scheme in this limit, as in the standard DQM. 
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Now, exploiting the gauge invariant first-class effective Hamiltonian, we perform the 
Legendre transformation to integrate out all the involved momenta, via the partition function 

r 3 2 _ 

Z = I Y[d qi d Pi dedne JJ S{Cl a )S(T b )det \ {Q a ,T b } | e*/* L , (5.5) 

i=l a, 6=1 

where T a (a = 1, 2) is a gauge-fixing function and the effective Lagrangian is given as 

L = Pi qi + ireO - H'. (5.6) 

Integrating out the momenta ttq and pi and exploiting the corresponding equations of motion 
and the constraints Q a , we can obtain the desired first-class Lagrangian 

L = L + L wz , 

Lwz = (2gi - q 2 ~ 26 + 3)6 - U\ (5.7) 

where Lq is given by Eq. Q3.1p . This Lagrangian ( |5.7[) is invariant under the transformation 

5q 1 = -e 2 , 5^2 = ei, 5g 3 = ex + 2e 2 , 56> = -e 2 , (5.8) 

with e 1 = €\ + 4e 2 and e 2 = — e\ — 2e 2 , which is obtained from the definition of 5c = {c, Q} 
with the symmetry generator Q = e a Q a . On the other hand, one can easily see that the 
canonical momenta (|3.2|) are modified in this first-class system as 

p 1 = q l + 26, p 2 = |(g 3 - q 2 ) + q-y + q 3 - 6, 

P3 = |(<?2 - qs), TTg = -6- 26, 

to yield the modified nonholonomic constraint Q 2 m terms of the qi and 6, 

Q 2 = -q 1 + q 2 -q 3 -2q 3 -l + 6, (5.10) 
which reduces into Eq. (|3.11 ) in the vanishing limit of the auxiliary field 6. 



(5.9) 



B. BRST invariant NHCS 



In this section we will obtain the BRST invariant Lagrangian in the framework of the 



Batalin-Fradkin-Vilkovisky formalism [IT| which is applicable to theories with the first- 



class constraints by introducing two canonical sets of ghosts and anti-ghosts together with 
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auxiliary fields (C a , V a ), (V a ,C a ), (N a , B a ), (a = 1, 2) which satisfy the super-Poisson algebra 
{C a ,P b } = {V a ,C b } = {N a ,B b } = 5 6 °.Q 

In the NHCS, the nilpotent BRST charge Qb, the fermionic gauge fixing function \1/ and 
the BRST invariant minimal Hamiltonian H m are given by 

Q B = C a Q a + V a B a , V=C aX a + V a N a , H m = H'-C l P 2 (5.11) 

which satisfy the relations {QB,H m } = 0, Q% = {Qb,Qb} = 0, {{^,Qb},Qb} = 0. The 
effective quantum Lagrangian is then given with H tot = H m — {Qb, ^} as follows 

L eff = P iq t + n 9 + B a N a + P a C a + C a V a - H tat . (5.12) 

Now we choose the unitary gauge x 1 = X 2 = ^2 and perform the path integration 
over the fields B u N 1 , C x , P 1 , Pi and C 1 to yield 

L ef f = Pi% + n e 9 + B 2 N 2 + P 2 C 2 + C 2 P 2 + P 2 V 2 - 7i e N 2 

-\{P\ - 2#D -qi-Q2- ql - (-4p 3 + 4 ?3 - 1)0 + ^tt 2 

+ (2p 3 -Pi- 2q 3 - 1)(N 2 + B 2 ). (5.13) 

Next, using the variations with respect to Pi, irg, V and P, one obtain the relations 

pi = q x - N 2 - B 2 , p 2 = ±(g 3 - g 2 ) + Ql + q 3 + 6> + TV 2 + B 2 , 

P3 = |(?2 - g 3 ) - 29 - N 2 - B 2 , ne = -9- iV 2 , (5.14) 

V 2 = -C 2 , V2 = t 2 , 

to, with the choice of iV 2 = — B 2 + 29, yield the effective Lagrangian 

Leff = Lq + L wz + L gh , 

L gh = - 1 -{B 2 ) 2 - {29 + 9)B 2 + C 2 C 2 , (5.15) 

^ere the super-Poisson bracket is defined as {A, B} = fylrffli - (-^) VAVB ^ \rj£\l where r] A 
denotes the ghost number in A and the subscript r and I the right and left derivatives. 
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which is invariant under the BRST transformation 

<W = -AC 2 , S B q 2 = AC 1 , <5 B g 3 = A^ 1 + 2C 2 ), 
M = -AC 2 , 6 B C a = -XB a , 5 B C a = 6 B B a = 0, 



(5.16) 



with C 1 = C 1 + 4C 2 and C 2 = — C 1 — 2C 2 , which are generalized transformation rules of Eq. 
( p.8| ), including the ghost fields. Here one notes that the first-class nonholonomic constraint 
Q 2 in Eq. ( |5.10|) is now generalized to include the ghost term contributions as follows 



^2 = -qi + Q2 - ?3 - 2 <?3 - 1 + 9 + B 2 



(5.17) 



and fli in Eq. (" |5. 1|) is trivially satisfied via using the relations (|5.14j) 



VI. CONCLUSION 

In conclusion, using the Hamilton- Jacobi (HJ) quantization scheme, we have investigated 
the nonholonomic constrained system, which possesses the structure of second-class con- 
straints, to compare with the standard and improved Dirac quantization methods (DQMs). 
We have shown that the integrability conditions in the HJ scheme are equivalent to the invo- 
lutive relations for the first-class constrained system in the improved DQM by constructing 
the first-class Hamiltonian and the corresponding effective Lagrangian in the framework of 
the improved DQM. Furthermore, with this effective Lagrangian, we have also constructed 
the BRST invariant nonholonomic constrained system. Through further investigation it is 
interesting to apply the improved HJ scheme to the constrained field systems as well as 
constrained point particle ones. 
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